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ABSTRACT 

Strong magnetization makes the disks surrounding young stellar objects rotate at rates that are 
too sub-Keplerian to enable the thermal launching of disk winds from their surfaces unless the rate 
of gas diffusion across field lines is dynamically fast. This under- appreciated implication of disk 
magnetization poses a considerable challenge for disk-wind theory. 

Subject headings: stars: pre-main sequence - accretion, accretion disks - ISM: jets and outflows 
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1. INTRODUCTION 

The pioneering study by Blandford & Payne (1982, 
hereafter BP82; see also Chan & Henriksen 1980) opened 
the door to a physical understanding of the highly colli- 
mated jets that emanate from magnetized disks accret- 
ing onto central gravitating masses. Pudritz & Norman 
(1983) made the first application to the disks in star for- 
mation, while Uchida & Shibata (1985) performed the 
first numerical simulations. Collectively, the disk-wind 
community has contributed to much of what we know 
about the magnetohydrodynamic (MHD) processes of jet 
launching, acceleration, and collimation. 

Invoking strong magnetization to aid in the magne- 
tohydrodynamics of the gas above the disk has serious 
consequences for the gas contained within the disk. War- 
die & Konigl (1993, hereafter WK93), Ferreira & Pel- 
letier (1995, FP95), Ferreira (1997, F97), Casse & Fer- 
reira (2000, CFOO), and Salmeron, Konigl, & Wardle 
(2007, SKW07) performed the primary studies of some 
of the effects. The result has been a shrinkage of the 
viable parameter space to launch disk winds. In fact, 
unless the gas diffuses across magnetic field lines at dy- 
namical rates, the sub-Keplerian rotation of the surfaces 
of realistic disks in young stellar objects (YSOs) may 
make magnetocentrifugally-driven disk winds a severe 
challenge if their surfaces are cold. 

In §2, we give the basic equations from the analysis 
of Shu et al. (2007, S07). In §3, we use this analy- 
sis to demonstrate that fast disk winds that are more 
than lightly loaded will not arise if they are thermally 
launched. In §4, we review how WK93, FP95, F97, 
CFOO, and SKW07 formally escape the conclusions of 
§3 by assuming that the diffusion of gas across magnetic 
field lines inside the disk occurs on dynamical time scales, 
and we discuss the astrophysical difficulties engendered 
by such a point of view. In §5, we summarize how X-wind 
theory surmounts the same difficulties. 

2. BASIC EQUATIONS 

Use a cylindrical coordinate system (m, tp, z) centered 
on a star with mass Af„, and consider centrifugal bal- 
ance in a thin disk rotating at radius vj with an angular 

Electronic address: fshu@physics.ucsd.edu 



velocity Q(w) that is a density-weighted average over z: 
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where G is the universal gravitational constant, £ is the 
surface density, and B z and B^ are, respectively, the 
midplane and upper surface values of the vertical and 
radial components of the magnetic field. From Ampere's 
law in cgs units, cB^ /2ir equals the vertically integrated 
current density J v , where c is the speed of light. Equa- 
tion |T]) then assumes that the Lorentz force per unit 
area, J v B z /c, associated with magnetic tension, provides 
the only opposition, apart from inertia, to the attraction 
of stellar gravity. We have ignored the self-gravity of the 
disk, which would additionally hinder the launching of 
magnctoccntrifugal winds. When integrating the diver- 
gence of the Maxwerll stress tensor over the disk thick- 
ness, we have also ignored terms like — (l/8Tr)d(B 2 )/dm 
as being a factor of the disk aspect ratio <C 1 smaller than 
the retained term {l/A-K)d{B z B w )/dz for a spatially thin 
disk. Retention of the effects of magnetic and gas pres- 
sure, which generally decrease with increasing radius w, 
would only worsen the problem of sub-Keplerianity de- 
tailed below. 

With / denoting the z-averaged fraction of Keplerian 
rotation at zu, ft = /T2k where fix = (GM^/w 3 ) 1 ' 2 , we 
obtain (S07, Appendix C) 
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Magnetic effects also affect the vertical structure of the 
disk. We define a characteristic thermal velocity a such 



that 



2zoP(0)/S, where P(0) is the midplane gas 



pressure. Integrating the equation of vertical hydrostatic 
equilibrium in z gives (S07, Appendix C): 
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where we define A and Aq as, respectively, the actual and 
thermal aspect ratios of the disk, 
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In the above, zq is a vertical height of the disk defined so 
that the z-integrated value of zp from the midplane to 
the surface of zero gas-pressure equals (zo/2)(£/2). Note 

that Aq = y/2(a/ tuflx) <C 1 in a thin disk where zuflx 
is the Kepleran speed at radius w. 

WK93, FP95, F97, CF00 r and SKW07 never employ 
the relation ^ to estimate /. Instead, they allow depar- 
tures from hydrostatic equilibrium above the disk that 
involve advective accelerations in the radial and verti- 
cal directions that are O(Aq) when measured in a non- 
dimensional sense. Section 4 demonstrates that such de- 
partures occur only if there is a transonic diffusion of gas 
across magnetic field lines that include a strong compo- 
nent of B v not present in the analysis of S07. 

Denoting B+/B z = I t with B z oc in self- 

similar models, and substituting equation §3§ into equa- 
tion ([2]), we get 



1 _p 4) Mo -1 



It \A A 



(5) 



The presence of a radial field £? ro that increases from 
zero in the midplane to B^ at the surface compresses 
the disk and makes A smaller than Aq. In the process, 
1 — f 2 becomes greater than zero. 

To make contact with disk-wind theory, let us follow 
FP95 in defining a "magnetization," 
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which also equals a quantity that WK93 call q?q. Dividing 
equation ([3]) by equation (O now yields I 2 = (2/p)(l — 
A 2 /Aq), where we have replaced B+/B z by It. Using 
this relation to eliminate A/Aq from equation §5§ now 
gives 
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where the positive quantity Ct(fi) = A/Aq is Ct(p) = 
(1 - /^/2)i/2 < i. 

The maximum value of p allowed mechanically in equa- 
tion ||7J) comes by setting / = 0. The aspect ratio then 
assumes the minimum value: 
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A mech /A , we obtain p mcch : 
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For p > /i m ech, magnetic forces are too strong for cen- 
trifugal equilibrium to be possible. For p = p m ech in thin 
disks where Ao/Ig <C 1, the roots |(5J| and ^ take the 
approximate forms, A mech « Aq/Ii, p mcc h ~ 2/ If. 

In the model of BP82 and in cases where the wind is 
relatively lightly loaded (e.g., F97), I = 1/4 or nearly 
so, i.e., B z oc nj~ 5 / 4 . The magnetic configuration just 
above the disk is then well approximated by a vacuum 
field, where the analysis of S07 yields It = 1.43. The 
corresponding inclination angle of the field from the ver- 
tical is i = 55°. With I t = 1.43, /x moch « 0.978. For p 
not closely approaching ^ mcc h, say, p sa 0.731 where the 
magnetic compression C(p) ~ 0.5, the right-hand side 



of equation |7]| is O(A ) <C 1. The analogous feature in 
numerical calculations, combined with the ameliorating 
effects of rapid diffusion, led disk- wind theorists to as- 
sert that the departure from Keplerian rotation remains 
small enough for other effects to overcome the deficit. 

3. THERMAL LAUNCHING 

Thermal pressure cannot do the requisite job in the 
absence of fast (turbulent) diffusion (see §4). Published 
models have the deficit of / from unity at the compressed 
surface z — z both larger than average (e.g., Fig. 5 of 
FP95) and smaller than average (e.g., Fig. 1 of SKW07). 
As a fiducial case, we assume the disk's surface to have 
the z-averaged / obtained from equation ([7]). To have a 
significant pressure exhaust of gas onto open field lines, 
we then require the thermal speed squared a 2 , at the sur- 
face of the disk to be (S07, eq. 31): 
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where we applied equations ((4| and ([7|) in the last step. 

Without external heating, the surface temperature of 
the disk, ~ 2~ x / 4 times the effective temperature, will 
be lower than its characteristic interior value, i.e., a 2 — 
Q s a 2 where O s is a fraction < 1. Equation (flQ|) then 
yields a quadratic equation for Cf (/zthcrm) = ^therm/Ah 
whose solution is 
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From CV(^therm) = AthermMo, we obtain ^thc 
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For thin disks where Ii<d s Ao -C 1, equations (fTTj) and 
(O give Ahcrm ~ A , /i t hcrm ~ 46 s A //f. The right- 
hand side of the second relation is, at most, several per- 
cent in the inner tenth of an AU in YSO disks. The 
magnetization is then so weak that the associated depar- 
ture from Keplerian rotation from equation ([7]), 1 — f 2 
ps 28 s j4oj would be similar to that implied by the radial 
gradients of the thermal gas pressure, i.e., very small. 
But weak magnetization, /i < /ithcrm, with good field 
loading, does not provide a long enough lever arm to 
launch fast disk winds even if / = 1 (see §3.5 of CF00). 

If the disk magnetization fi is chosen greater than 
Mtherm, but less than /imoch, the actual departure 1 — f 2 
~ O(Aq) exceeds the depth of the dimensionless effective 
potential 1 — f 2 ~ O(Aq) that the surface temperature 
can offset. Because a 2 may be less than 1% of GM^/w 
in the inner disk, even / = 0.99 may constitute a sub- 
stantial barrier by these standards. As argued earlier, 
no room exists for a compromise p between /Xtherm and 
Mmcch that will drive disk winds that both load and fling 
if one depends on thermal pressure to launch. 

4. DIFFUSIVE LOADING 

WK93, FP95, F97, CF00, and SKW07 differ with the 
analysis of §§2 and 3 because they do not invoke thermal 
pressure to launch disk winds. Instead, they overcome 
local deficits in / by diffusive loading onto magnetic-field 
lines coupled to magnetocentrifugal acceleration. 
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For an imperfectly conducting, lightly ionized, gas, suf- 
ficiently collisional to act as a single-component fluid of 
velocity v, the induction equation reads 



dB 

~dt 



V x (B x v) = —V x [rj ■ (V x B)] , 



(13) 



where 77 is a generalized tensor resistivity associated with 
Ohm's law (Norman & Heyvaerts 1985, see especially 
eq. 22). In a frame that corotates locally with the mag- 
netic held, so that dH/dt = 0, we may "uncurl" equation 
(fl3|) and obtain the balance of advection and diffusion, 
B x u = —rj- (V x B), where u is the fluid velocity relative 
to the rotating held lines. 

Keeping only the dominant terms in a thin disk, we 
follow WK93 and SKW07 in decomposing the advection- 
diffusion balance into its tp and w components: 
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where the dimensionless rotation of the held fg is a func- 
tion of z at fixed zu (in such an Eulerian description, we 
are not following a given held line). For WK93, r\\ is as- 
sociated with the coefficient of ambipolar diffusion, and 
772, with the inverse of the Hall conductivity. In F97, and 
CF00, 771 is replaced by v m in equation (fl4|) and by 
in equation (| 15(1 . with the other, non-diagonal, coefficient 
772 set to zero. 

In the disk proper (where u z is small) of both Figure 
1 of SKW07 and Figure 5 of FP95, the terms on the 
right-hand sides of equations (fl4|) and (|15|) contribute, 
respectively, negatively to the radial drift and positively 
to the tangential drift relative to the held lines. These 
relative motions reverse as the gas climbs vertically into 
the disk wind. Eventually, the right-hand sides vanish 
as ideal MHD is approached (771 and 772 — > 0), with u 
becoming parallel to B in the corotating frame. The 
physical interpretation of the result is as follows. 

A held line rotates at a dimensionless angular veloc- 
ity /b that is intermediate between the rotation rate / 
of the gas near the midplane layers of the disk and the 
gas in the upper atmosphere that blows into a disk wind. 
This rate is fixed in steady state so that angular mo- 
mentum balance is reached following a given held line 
(a condition not examined in this paper). In the WK93 
description, the neutral gas near the midplane would ro- 
tate at the Kepler speed except for collisions with ions 
and electrons that are tied to the more slowly rotating 
magnetic held lines. The collisions therefore make the 
neutral gas lose angular angular momentum, which then 
cause the neutrals to drift radially inward relative to the 
radially stationary held lines and the ions and electrons 
tied to them. 

Freezing to held lines in the upper atmosphere of the 
disk yields the familiar effect of a "sliding bead on a 
rotating rigid wire" (Chan & Henriksen 1980) that ac- 
celerates the gas through the slow MHD transition as 
it gains angular momentum at the expense of the mat- 
ter in the deeper layers of the disk. Thus, the gas that 
participates in the disk wind ends up rotating at super- 
Keplerian speeds. A similar phenomenon occurs in the 



FP95 description, except that the role of collisions be- 
tween neutrals and ions/electrons is replaced by a mixing 
of fluid parcels created by MHD turbulence. 

When an outwardly bending field line rotates uni- 
formly along its length, there will always be a point where 
centrifugal effects on an electrically conducting test par- 
ticle placed there counterbalances an opposing gravita- 
tional field. In Figure 5 of FP95 and Figure 1 of SKW07, 
this point occurs at ~ 2vt/^k, where vt is the thermal 
velocity in the midplane. For a vertically isothermal disk 
(the most optimistic case), where a — (2/y / 7r)wT, this 
point is ~ 2.5zq since zo ~ 0.5Aow ~ 0.8vt/^k for typ- 
ical parameters used to launch disk winds. The height z$ 
is dchned by where the condition of vertical hydrostatic 
equilibrium would have created a surface. As in any at- 
mosphere, the gas density does not really drop abruptly 
to zero at z$, but unless diffusion works very efficiently, 
there cannot be much gas at 2.5z . 

Trouble on this point enters when numerical solutions 
are adjusted to give large rates of diffusion. In order 
of magnitude, let us estimate dB^/dz ~ B+/zq and 
dB v /dz ~ B+/z in the disk proper. In the more suc- 
cessful cases considered by FP95, the quantities 771 ~ 
z a ~ 772- With £?+ > and — B+ > comparable in 
magnitude to B z > 0, equations (TTJJ and (Ti"5)) , with u z 
negligible in the disk proper, now yield the estimates, 

u ro a, f - f B Q ~ A . (16) 

y CrM*/m 

The first deduction, that the inward drift speed in the 
bulk of the disk approaches sonic values, is verified, for 
example, by Figure 4 of FP95 where u m » — vt- The 
second deduction, that the dimensionless rotation rate 
of the gas / can exceed that of the held fg by O(A ), 
would produce favorable conditions for the launching of 
disk winds if the field rotates slower than the gas in the 
disk proper, but still within fraction O(Aq) of Keplerian. 

In the work of FP95, the diffusion coefficients 771 and 772 
are chosen large enough so that, unlike S07, the magnetic 
field is not well coupled to the matter in the deep interior 
of the disk. With 771 ~ zqci, the diffusive time scale to 
straighten poloidal held lines that are bent on a scale zq 
is ~ Zq/t]i r~j Zo/a, which is of order the inverse Kepler 
rate fi^ 1 (or shorter, if the disk is severely magnetically 
compressed). The dynamical slippage of the matter past 
the field, both in the w and ip directions, is a serious 
weakness of the models. 

Rapid diffusion can carry gas to the slow MHD transi- 
tion point beyond which magnetocentrifugal fling works 
in the outward direction to produce an outflowing wind. 
Such a strategy generates, through equations (|14[) and 
(Ti7j|) , radial drift speeds inside the disk that are transonic 
~ — vt- Consider what this implies if a disk wind were 
responsbile for the primary accretion mechanism of the 
disk all the way out to an outer disk radius ~ 300 AU in 
a low-mass YSO. The midplane temperature at 300 AU 
is unlikely to be much less than, say, 20 K, which corre- 
sponds to an isothermal sound speed in cosmic molecular 
gas of vt = 0.278 km s™ 1 . The whole disk would then 
shrink onto the star in 300 AU/(0.278 km s" 1 ) - 5,000 
yr, which is far too short a time scale to be realistic. 

If global magnetorotational instability is a viable al- 
ternative for YSO accretion disks, the analysis of S07 
and Shu et al. (2008) implies that the turbulent resis- 



4 



tivity rji is smaller than estimated by FP95 or F97 by 
an additional factor ~ [iA <C 1. If we continue to as- 
sume that — B+ and B+ are comparable to B z , with the 
former holding only if disk winds can be launched, then 
equations (fl"4|) and (| 1 5[) imply 

u^^-fiAa, f - f B ~ fiAA < A%. (17) 

The second relation above states that diffusion can now 
make up deficits of fg only by amounts < Aq, where 
/b itself is given by / to O(Aq) in such circumstances. 
In turn, / departs from Keplerian rotation typically by 
O(Aq) for strongly magnetized disks according to equa- 
tion ([7]). A low (S07; see also Lubow et al. 1993) rather 
than a high (FP95 or F97) rate of field diffusion yields 
equation (fT0|) as, after all, the correct approximate cri- 
terion for wind loading and launching. 

The numerical choices made by WK93 are more con- 
servative than those of FP95, and the net results for 
u ro and / — fg are intermediate between the high- and 
low-diffusion scenarios described above (see Appendix 
of WK93). Specifically, in the numerical case depicted 
in Figure 1 of SKW07, where ao = ^/Jl = 0.95 and 
a 2 s = v% = Q.QIGA'h/vj, the dimensionless loading is 
k = 3.2 x 1CP 4 and the square of the lever arm is A = 395 
(these are and J w , respectively, in the notation of X- 
wind theory), corresponding indeed to a disk wind that 
can fling, but is very lightly loaded. Such winds do not 
describe observed YSO jets well. The order of magnitude 
changes (exponential sensitivity) in quantities like k and 
A when rj changes by a factor of a few gives powerful tes- 
timony to the crucial role that magnetic diffusivity plays 
in current theories of the launching of disk winds. 

5. DISCUSSION AND CONCLUSION 

Disk winds have been invoked as a solution to the angu- 
lar momentum problem in accretion disks (e.g., Pcllcticr 
& Pudritz 1992). The challenge of sub-Keplerian rota- 
tion facing magnetocentrifugally-driven, cold, disk winds 
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for much, if not most, of its radial span suggests that in- 
ternal mechanisms other than the back reaction to disk 
winds must account for the primary accretion mechanism 
of such disks. Examples of such mechanisms include spi- 
ral density waves for high-mass disks and the magnetoro- 
tational instability for low-mass disks. 

Disruption of the disk at an inner edge by a stellar 
magnetosphere will drive an X-wind (Cai et al. 2008 and 
references therein) if the resistivity in the disk is appre- 
ciably smaller than the viscosity (for numerical simula- 
tions, see Romanova et al. 2008). How does X-wind the- 
ory escape the conundrum posed by equation |T])? Near 
the inner edges of disks, the X field geometry involves a 
swing of outwardly bending field lines (right half of X) to 
inwardly bending field lines (left half of X) as accretion 
disks are truncated by funnel flows through their interac- 
tions with a stellar magnetosphere. The swing of B+B z 
from positive values through zero (vertical field lines that 
are "dead" to outflow or inflow) automatically promotes 
a transition from sub-Keplerian to Keplerian rotation. 
Indeed, when the density decreases sharply inward on 
a fractional radial scale ~ Ao, the no-longer ignorable, 
extra, radial push of gas pressure toward the inner edge 
speeds up disk rotation and makes the footpoints of some 
of the outwardly bending field lines exterior to the dead 
zone sufficiently close to Keplerian rotation as to enable 
the thermal launching of X-winds (see the discussion of 
Shu et al. 2008). 
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